Conductance of interacting Aharonov-Bohm systems 
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A simple formula for the zero-temperature linear response conductance of an interacting meso- 
scopic region, threaded by magnetic flux, and attached to noninteracting single-channel leads is 
presented. The formula is valid for a general interacting system exhibiting Fermi liquid proper- 
ties. As an example of the efficiency of the formula the results for the conductance of a simple 
Aharonov-Bohm ring with Kondo-Fano resonance physics are presented and compared with numer- 
ical renormalization group results. 



The electron-electron interaction often plays a cru- 
cial role in transport through mesoscopic systems. The 
Kondo effect in quantum dots is e.g. a prototype 
phenomenon, where correlations dominate the conduc- 
tancei* 2 -. Another remarkable effect is the Fano reso- 
nance physics intertwined with Kondo physics^. Recent 
advances in nanofabrication of Aharonov-Bohm (AB) de- 
vices made possible a realization and a systematic study 
of such phenomenal. 

For systems where the interaction is absent, or very 
weak, transport properties may be determined using the 
Landauer-Biittiker formalism^. On the other hand, if the 
interaction is important in the system being studied, the 
use of a more general approach is essential. An appro- 
priate formalism, expressing the conductance in terms 
of non-equilibrium Green's functions, was developed by 
Meir and Wingreeri The formalism can in principle 
be used to treat systems at a finite temperature, finite 
source-drain voltage, and can even be extended to de- 
scribe time-dependent transport phenomena^. Another 
approach to describing the transport in interacting sys- 
tems, applicable only to systems in the linear-response 
regime, is the Kubo formalism^. 

Recently, a simple method for calculating the conduc- 
tance through a region with electron-electron interaction 
(e.g., a molecule, a quantum dot, a quantum dot array 
or a similar 'artificial molecule' system, ...), connected 
to single-channel noninteracting leads, was presented in 
Ref.licl (hereafter referred to as RR)^i. The conductance 
is determined solely from the ground-state energy of an 
auxiliary system, formed by connecting the ends of the 
leads of the original system into a ring and threaded by 
magnetic flux. The method is applicable to Fermi liquid 
systems at zero temperature and in the linear response 
regime. The validity is additionally restricted to the 
class of interacting systems obeying time reversal sym- 
metry. Several conductance formulae, becoming exact 
in the limit of a very large ring, are derived, the most 
practical being the two-point formula, 
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where E (0) and E (n) are the ground-state energies of 
the auxiliary interacting system with periodic and an- 



Figure 1: (a) Interacting mesoscopic region (gray-shaded 
sites), threaded by magnetic flux and coupled to noninteract- 
ing leads, (b) Auxiliary ring system, (c) Behavior of energy 
levels as the flux threading the ring is varied, (d) System from 
Ref.H with a quantum dot embedded in an AB ring. 



tiperiodic boundary conditions, respectively. The aver- 
age level spacing at the Fermi energy A = [Np (e^)] 1 is 
determined by the density of states in an infinite lead p (e) 
and the circumference of the ring TV, and Go = — is the 
conductance quantum. The corresponding ground-state 
wave function can be chosen real, which simplifies the 
determination of the ground-state energy by, for exam- 
ple, variational or density matrix renormalization group 
approaches. 

In this paper we present a generalization of the conduc- 
tance formula Eq. to systems which explicitly exhibit 
time reversal asymmetry, such as is e.g. an AB-type of 
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system presented in Fig.nja). The formalism presented 
here is valid for general systems consisting of intercon- 
nected sites with interaction, threaded with magnetic flux 
and connected to noninteracting leads. 

First we consider the noninteracting case and single- 
electron states. The transmission amplitudes t k and t' k 
of an electron with a wave vector k and an energy e k de- 
scribe the transmission from the left to the right lead and 
from the right to the left lead of the system presented in 
Fig. da), respectively. They are related by an expression 
derived in RR, 
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where $ = §0 is the magnetic flux through the auxiliary 
system shown in Fig.njb). If there is no AB flux thread- 
ing the mesoscopic region, the time reversal symmetry 
is restored, tk = t' k and the energy is an even function 
of (f>, as illustrated in Fig. ^c). In a general case, the 
unitarity of the scattering matrix requires \tk\ = |tl| and 
the transmission amplitudes are related as tk = t k e l ^ ok 
and t' k — t k e~ l ^ ok . Expressing tk in terms of its modulus 



\tk\ and a phase-shift (p k , tk 
equation Eq. <|2jl reads 



\tk\e ltfk , the eigenenergy 
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Algebraic manipulation equivalent to that in RR gives an 
implicit equation for |i(efc)|, exact to the leading order 
in jj, 

da,rccos(T\t(ek)\co$[(f)-(f>o(£k)]) ,, , N de k 
= TrNp{e k )- 
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where the sign =F depends on weather k belongs to a 
decreasing (+) or to increasing (— ) branch of the cosine 
function in Eq. 10 and p{e) = {ir^jAt^ — e 2 ) with to 
being the nearest neighbor hopping in the leads. By using 
the principle of mathematical induction, as in RR but 
generalized to a finite (j>o , the transmission probability at 
the Fermi energy Ef is expressed as 
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where E is the ground-state energy of a system contain- 
ing an even number of electrons. Neglecting the vari- 
ation of <(e_f) and p(sf) with 4> (the error made is of 
the order of -^), Eq. iJ^J leads to a universal form of 
the ground-state energy of the auxiliary system E((f>) 
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this form, the transmission probability can be extracted, 
and the conductance is given by 
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where 4>o = (f>o ( £ f) is determined by the position of the 
minimum (or maximum) in the energy vs. flux 4> curve, 



schematically shown in Fig. ^c). The conductance can 
also be calculated from the more convenient four-point 
formula^: 
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and 4>o can be determined from the expression 
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We now consider the interacting case. If no AB flux is 
present in the mesoscopic region and the ground state of 
the system exhibits Fermi liquid behavior, i.e., the per- 
turbation theory in the interaction strength is valid, the 
imaginary part of the self-energy due to the interaction 
vanishes quadratically at the Fermi energy, 



ImEjj (cj + id) oc (uj — e f) , 



(9) 



and a quasiparticle Hamiltonian can be constructed for 
each value of ef- As shown in RR the conductance cal- 
culated from the quasiparticle Hamiltonian reproduces 
the conductance of the interacting system, and is exactly 
given by Eq. [or Eq. since E (tt/2) = E (-7r/2) in 
this case]. If the time reversal symmetry is broken due 
to AB flux, Eq. © is not valid and the proof has to be 
reconsidered. Repeating the steps as presented in detail 
in RR, the proof is restored and basically unchanged if 
the exact self energy obeys the relation 



t [Ejj (u> + id) — Ey (l) — iS)\ oc (u> — ep) 2 ■ (10) 



It follows that also the linear response conductance of an 
interacting AB system at zero-temperature is given by 
the four-point formula Eq. {JJ . This condition is fulfilled 
if the system is a Fermi liquid. 

In order to demonstrate the practical value of the 
method, we quantitatively analyze the conductance 
through an Aharonov-Bohm ring with a quantum dot 
placed in one of the arms as presented in Fig. njd) 3 ' 4 . 
The quantum dot is described as an Anderson impurity 
with level position ed and a charging energy U, coupled 
to each of the leads with a tunneling matrix element t\ . 
Electrons can also be transferred from one lead to the 
other directly through the upper arm of the AB ring. 
This process is described by a tunneling matrix element 
ti. The ring is threaded by an AB flux <&ab = \§ab in 
such a way that only the direct tunneling matrix element 
is affected, i.e. ti — » t2e lc ^ AB . The Fermi energy is set 
at the middle of the band, thus A = ^y 2 -. We perform 
the finite-size analysis of the four-point formula Eq. JJJ , 
changing the circumference of the ring N . In order to 
be able to compare our results with those of the NRG 
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Figure 2: The conductance of a noninteracting (U = 0) AB 
ring at 4>ab = f • The inset shows convergence of the con- 
ductance at the peak near Ed = 0. Parameters: ti = 0.177to, 
h = 0.298to. 



method, we choose the same values of parameters as in 
Ref.0. 

In Fig. the results for the conductance of a noninter- 
acting (U — 0) system with 4>ab = f are presented. Due 
to a non-zero AB flux, the ground-state energy must be 
determined at four <fi points as required by Eq. 0. The 
conductance exhibits a typical Fano resonance with a dip 
and a sharp peak. Results calculated for various numbers 
of sites in the ring N are compared with the exact con- 
ductance curve. The inset shows a convergence test of 
the method at the peak near td = 0, where due to a 
strong energy dependence of the transmission amplitude 
the convergence is the most delicate. 

In the interacting case, we determined the required 
ground-state energies using an approach similar to the 
variational method of Gunnarson and Sch6nhammei»ifl»A^ 
for the Anderson model. The variational basis set is gen- 
erated from the ground state 1 0} of an auxiliary noninter- 
acting Hamiltonian in which the energy level of the dot Ed 
and the hopping matrix element t\ between the dot and 
the leads are renormalized in such a way to minimize 
the ground-state energy. Apart from the Hartree-Fock 
case where this is the only variational wave function, we 
form two additional basis sets. The first has three varia- 
tional wave functions Pq |S), P\ |6) and Pi |5), where Pq, 
Pi and P2 are projectors onto unoccupied, singly occu- 
pied and doubly occupied dot site. In the second basis 
set we add four additional wave functions, PqVPi |6), 
P 2 VP\ |0), PiVPq j 6) and P1VP2 |6), where V is the op- 
erator describing hopping between the leads and the dot. 
In all the three cases, the noninteracting limit is cor- 
rectly reproduced. Furthermore, except in the Hartree- 
Fock case, the method is also exact in the limit where 
the dot is decoupled from the rest of the system. 

An interacting AB ring with U ^ serves as a good 
nontrivial test of the method. We choose a strong cou- 
pling regime as in Ref. H with U = 8r, where T = 
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Figure 3: Zero-temperature linear response conductance as a 
function of level position Ed for various values of the AB flux. 
Labels '3' and '7' correspond to variational wave functions 
with three and seven terms, respectively. The dash-dotted 
line is the result obtained using the HF approximation and 
the dotted line is the NRG result from Ref. Q, Parameters: 
h = 0.177to (r = 0.125t ), ti = 0.298t , U = t = 8T. 



47rffp (ef) is the line width of the dot level in the ab- 
sence of the upper arm of the AB ring. The results are 
presented in Fig. and compared with results obtained 
using the NRG method^, together with the correspond- 
ing Hartree-Fock (HF) curve. It should be noted that 
with the parameter set chosen, the underlying physics is 
in the strong correlation regime, where e.g. the width of 
the Kondo peak is of the order 10~ 4 U. The HF method 
therefore fails to reproduce NRG results, while the vari- 
ational results are very close to the NRG curve. By 
extending the variational space from 3 to 7 terms, the 
agreement with NRG improves even further and the small 
remaining discrepancy would most probably be addition- 
ally reduced if an even richer variational wave function is 
used. Interestingly, we find the largest deviation in the 
'empty orbital' regime, Ed < —U. The results of the vari- 
ational and the HF methods here agree, both deviating 
from the NRG result, while asymptotically all the curves 
reach the correct limit G = 0.3G(r». We have also checked 
the result in this regime using the second order pertur- 
bation theory approach^ which gives a conductance in 
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Figure 4: The same as in Fig.^b), but calculated for different 
numbers of sites in the auxiliary ring N. The inset shows 
convergence of the conductance at td = —(7/2. 



agreement with variational and HF methods (not shown) . 

In Fig.Q]a convergence test of the method for the result 
from Fig. Efb) is shown. The convergence with N is fast 
in the empty orbital regime and becomes progressively 
slower as td shifts toward the Kondo regime. The reason 
can clearly be attributed to a very strong coupling regime 
resulting in an extremely narrow Kondo peak. A very fine 
energy resolution is therefore required to resolve the peak 
and to obtain a converged conductance curve. 

Broken time reversal symmetry in AB systems is sig- 
naled by 0o 7^ 0. In Fig. is shown the phase shift <po 
corresponding to Figs. El and Efb) as determined from 
Eq. (jSJ. In contrast to the smooth noninteracting result, 
the phase shift in the interacting case exhibits a well de- 
veloped plateau, corresponding to the Fano-suppressed 
Kondo plateau in the conductance. The (j)o curve, which 
clearly cannot be correctly reproduced in the HF approx- 
imation, can be explained as follows: In the empty or- 
bital regime, the current mainly flows through the upper 
arm of the ring and therefore, electrons acquire an addi- 
tional phase shift 4>q ~ 4>ab (note that <fio and </>o — 7r are 
physically equivalent). On the other hand, in the Kondo 
regime almost all the current passes through the quan- 
tum dot and no additional phase shift is present, 4>q ~ 0. 

In summary, we have derived a formula for zero- 
temperature linear response conductance of an interact- 
ing electron region coupled to single-channel leads where 



the electron-electron interaction is absent. The interact- 
ing part of the system can be a general Aharonov-Bohm 
type of interferometer with broken time reversal symme- 
try [Fig, HJa)]. The conductance of such an open system 
is exactly determined from an auxiliary closed system 
[Fig-IHb)], where the leads of the original system are con- 
nected to form a ring of N noninteracting sites threaded 
by an auxiliary magnetic flux. Eq. Q, which follows from 
the universal form of the ground-state energy of the aux- 
iliary system for a large but finite N, expresses the con- 
ductance in terms of the ground-state energy evaluated 
at four different values of the auxiliary flux. The proof 
of validity of the formula for interacting systems relies 
on the mapping of the system onto an effective quasipar- 
ticle problem and is therefore valid for systems exhibit- 
ing Fermi liquid properties. We have demonstrated the 
usefulness of the formula by applying it to a prototype 
system exhibiting Kondo-Fano behavior. Results based 
on the four-point formula and variational ground-state 
energies confirm results of the numerical renormalization 
group method. 
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Figure 5: Phase shift 4>o as a function of Ed /to for the non- 
interacting AB ring of Fig. |5] (NI) and for various ground- 
state energy methods applied to the interacting AB ring of 
Fig. lib). " 
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